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Abstract 

We have extended the high-temperature susceptibility series of the 

three-dimensional spin-i Ising model to 0(t> 26 ). Analysis of the new 

series gives a = 0.101 ± 0.004. 

O 



In an earlier paper || we gave series to order v 22 for the high-temperature ex- 
pansion of the zero-field partition function of the 3-dimensional Ising model. 



More precisely, we gave the coefficients a n , n = 0, 22, defined by 
Z = 2[cosh(J/kT)] 3 $(v), with $(v) = 5>„t/\ 

n 

The series were obtained by the finite-lattice method. One difficulty with 
the finite-lattice method for this problem is its voracious appetite for com- 
puter memory. Our earlier computation in fact calculated the series to two 
further terms - to order v 26 - but due to addressing limitations, we were 
unable to retain the intermediate information. This particular calculation 
requires 2.08GB of memory, and we were unable to address more than 2GB, 
due to operating system limitations. We have now been able to re-run our 



program under a different operating system that permits us to address this 
large address space. 

The program was run on an IBM 3090/400J with 500MB of memory 
and 2GB of extended storage - a slower type of memory. The use of the 
MVS operating system allowed the large address space to be used. Even so, 
2-byte integers were used, and the program run twice modulo two different 
primes. The results were combined using the Chinese Remainder Theorem, 
and provided the least significant digits of the new coefficients, the most 
significant digits were obtained by differential approximants. The final results 
were then compared by running with a third prime. Each run took 150 hours. 

As a result, we have obtained two further non-zero terms (the partition 
function being an even function has vanishing odd-order coefficients). We 
have also obtained the 6 most significant digits of the 0(v 28 ) coefficient, by 
the method of differential approximants. In our earlier paper we obtained 
the coefficient of the 0(v 24 ) coefficient by this method, and claimed the co- 
efficient to be a24 = 27337* 10 7 . The present calculation gives the coefficient 
as a24 = 273374177222, verifying our prediction. The subsequent coefficients 
are found to be a 26 = 4539862959852 and a 28 = 7474452 * 10 7 ± 5 * 10 7 , 
where the last coefficient is obtained by differential approximants. (The ap- 
proximate coefficient was not used in the subsequent analysis, as differential 
approximants require more accurate coefficients. It is nevertheless useful for 
ratio type methods of analysis). 

As we were completing this work we received a preprint Q in which a 
variant of the finite-lattice method, using helical boundary conditions was 
used to obtain one further coefficient than we had previously obtained. This 
work also confirmed our predicted coefficient, and agrees with our exact co- 
efficient. (Note that they give the free-energy series and we give the partition 
function series). They also predicted a 2 6, and our exact coefficient confirms 
their predicted value. 

We have analysed the new series by several methods. The series is now, 
for the first time, sufficiently long that the method of differential approxi- 
mants can be used with some confidence. For our initial analysis, we used 
unbiased approximants, but for maximum precision we used biased approxi- 
mants. This requires a knowledge of the critical temperature which has been 
accurately estimated from the more readily analysed high-temperature sus- 
ceptibility series, as well as from a variety of Monte-Carlo estimates. The 
series estimates are reviewed in H, and we use the best estimate given 



there, v c = 0.218093, which is in good agreement with the most recent, 
high-precision Monte-Carlo estimate, of v c = 0.2180992 ± 0.0000026. 

Our method of analysis is fully described in || , and provides a weighted 
mean of critical exponent estimates from inhomogeneous first- and second- 
order differential approximants, with one estimate obtained for each order of 
the series. Our analysis was carried out on the coefficients of the partition 
function itself. Our unbiased estimates are 

v 2 c = 0.04756 ± 0.00003 and 2 - a = 1.905 ± 0.016 with K = 1 

v 2 c = 0.04756 ± 0.00002 and 2 - a = 1.897 ± 0.012 with K = 2. 

In the above, K = 1, 2 refers to first- and second-order differential approxi- 
mants respectively. The unbiased estimates are seen to be in excellent agree- 
ment with the sussceptibility series estimate v 2 = 0.0475646, while an esti- 
mate of a = 0.10 ± 0.01 can be made. A biased analysis yields the following: 

2 - a = 1.899 ± 0.004, K = 1 and 2 - a = 1.900 ± 0.006, K = 2 

Thus we find, from this analysis, a = 0.101 ± 0.004. This is substantially 
more precise than our earlier analysis, using two fewer series coefficients, 
of a = 0.104 ± 0.018. It is consistent with the analysis of who find 
a = 0.104 ± 0.004, though as can be seen we favour a rather lower value. 
Note that second-order differential approximants implicitly take correction- 
to-scaling terms into account. The agreement between first- and second-order 
differential approximants suggests that correction-to-scaling exponents are 
weak. A subsequent analysis provides numerical confirmation of this. 

Ratio techniques can also be used with this series. We have analysed 
the free-energy series by a variety of extrapolation methods, based on the 
observation that if the free-energy, \1/ jkT ~ A{\ — v 2 /v 2 ) 2 ~ a , then the ratio 
of successive coefficients in the series expansion of ty/kT behaves like -4(1 + 

9 —^-), with higher order corrections from correction-to-scaling exponents, as 
well as corrections due to analytic terms. In any event, the sequence of 
ratios can obviously be re-arranged to give a sequence that will converge to 
a. Neville extrapolation (which takes into account only analytic correction 
terms), gives a = 0.103 ± 0.006. Other extrapolation methods, such as 
Levin's u-transform and Brezinski's ^-algorithm are less accurate, allowing 
only the estimate a = 0.10±0.03. 



In our previous analysis, we also studied the amplitude of the "correction- 
to-scaling" term, a$, where the specific heat is defined to have the scaling 

form C ~ A|t|~ Q [l + a e \t\ e + a x \t\ -\ ], where t = (T - T c )/T c and 9 « 0.52 

[§]. In [0] it was argued that ag should be negative, and our earlier analysis 
11 seemed to confirm this, in that we found a# ~ —0.04. This can be seen 
from the behaviour of the ratios of successive coefficients, as follows: We 
first write C(v) = X)c n t> 2n , as the expansion we obtain is in terms of the 
usual high-temperature variable v =tanh( J/kT). Note that, to leading order, 
t = (T — T c )/T c = B(v — v c )/v c , where B is a positive constant. It therefore 
follows that the correction-to-scaling amplitude of the specific heat series 
expanded in the variable v 2 should also be of negative sign. Writing 

C(v) = £ c n v 2n = A(l - v 2 /v 2 c r a (l + 6(1 - v 2 /v 2 c ) e + •■•), 

it follows that 

AT(a + n) . bT(a)T(a + n - 6) 

■H , \ ; H 



hence 



^ n T(a)T(n + l)v? n[ T(a-6)T(a + n) 






c n _i f 2 n r(a — 6)n 0+l n 2 

Taking a w 0.1 and ~ 0.5, it follows that the above equation can be 



rewritten as 



1 M a-1 1.28..6 _, 1 N1 



c n _] f 2 n r7, 6 ' +1 n 2 



Hence we find that 



c n 2 _ . 1.28. .6 „,1 



-v: - l)n + 1 ~ a + " + O(-). 

'c n -i n y n 

This means that if b < 0, estimators of a, given by the l.h.s. of the above 
equation, should approach a from below. In fact we find the approach to 
be from above, but a simple n-shift of 1 makes the approach change to an 
approach from below! Even an analysis taking into account the analytic 
correction term does not alter this behaviour. To be more precise, we have 
repeated the above analysis with an additional analytic correction-to-scaling 
term present, and found that the numerical value of b changes sign with an 
n-shift of just 1. In all cases, the estimate of b is numerically rather small, 



a estimate 


Method 


Reference 


0.101(4) 
0.104(4) 
0.1100(45) 
n 1 n4 + 0- 006 

U.1U4_ 0009 

0.110 


Series 

Series 

Field theory 

Series + hyperscaling 

Field theory + hyperscaling 


This work 

E 
E 

1 
E 



Table 1: Summary of a estimates 



and we conclude that this analysis is not sensitive enough to distinguish b 
from zero. A similar conclusion, based on a somewhat different analysis, was 
obtained in [p. 

Our estimate of a is rather lower than the field-theory estimate of a = 
0.110 ±0.0045, but the field-theory and series estimates are both (separately) 
consistent with the hyperscaling relation di> = 2 — a. Our best series estimate 
of v = 0.632±g;gg§ implies a = 0.104±g;gg|, while the best field-theory estimate 
U is v = 0.630, which implies a = 0.110, a value at the centre of the field- 
theory estimates. 

We summarise the various estimates of a in table p]. 
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